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Abstract 

We investigate the dynamical behavior of a non-ideal Duffing oscillator, a system com¬ 
posed of a mass-spring-pendulum driven by a DC motor with limited power supply. To 
identify new features on Duffing oscillator parameter space due to the limited power 
supply, we provide an extensive numerical characterization in the bi-parameter space by 
using Lyapunov exponents. Following this procedure, we identify remarkable new pe¬ 
riodic windows, the ones known as Arnold tongues and also shrimp-shaped structures. 
Such windows appear in highly organized distribution with typical self-similar structures 
for the shrimps, and, surprisingly, codimension-2 bifurcation as a point of accumulations 
for the tongues. 
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1. Introduction 


In recent years, there has been an increasing amount of work on nonlinear dy¬ 
namics characterizing the possible structures in two-dimensional control parameter (bi¬ 
parameter) space [lj. Accordingly, periodic windows with important features, mainly 
shrimp-shaped structures |2] and Arnold tongues (3;, 0,0, have been identified in several 
systems such as two-gene model i6|, impact oscillator [3, 


tivistic particles 0, tumor growth model 0 , Chua’s circuit 
model [10|, and Red Grouse population model [15|. 


dissipative model of rela- 


11J, 12), 13J, prey-predator 
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In the nonlinear dynamics context, oscillators with mechanical coupling have recently 
attracted a significant attention due to the complexity of the dynamics for high de gree - 
of-freedom devices and possible applications to advanced technologies [lfi |l7j, |l8l (19j, 20]. 
Among the class of mechanical coupling oscillators an interesting example is the mass¬ 
spring-pendulum system 2l|, [22|. Svoboda and collaborators studied a system of masses 
with a pendulum, where the pendulum is attached to one mass of a chain of masses 
connected by springs |0. They showed that auto parametric resonance can arise. In 
Reference 24] was investigated the influence of nonlinear spring on the auto parametric 
system. It was verified the existence of regular and chaotic motions. 

In this work, we investigate the parameter space organization of a non-ideal Duffing 
oscillator, namely, the mass-spring-pendulum system. Duffing oscillator is a forced oscil¬ 
lator with a nonlinear elasticity, and it is described by a non linear differential equation 
of second-order that has been used in a variety of physical processes. This oscillator is 
well known in engineering science, and it has been used to model the dynamics of types 
of electrical and mechanical systems. Almong and collaborators experimentally studied 
signal amplification in a nanomechanical Duffing resonator via stochastic resonance pfij ] . 
The Duffing oscillator is also a useful model to study the dynamics behavior of structural 
systems, such as columns, gyroscopes, and bridges [26]. 

The non-ideal character of the studied oscillator is a consequence of the fact that 
the source of energy is given by a DC motor with limited power supply 


27, 


. Previ¬ 
ous studies of this system have shown a rich dynamical behavior with several nonlinear 
phenomena, like quasi-periodic attractors, chaotic regimes, crises, coexistence of attrac¬ 
tors, and fractal basin boundaries 0, [3ol . Isij] . Here, our main purpose is to provide a 
global parameter analysis of the behavior of this oscillator with a mechanical coupling. 
The main features found in the parameter space were the self-similar structures, and 
codimension-2 bifurcation as a point of accumulations for the Arnold tongues. Compar¬ 
ing with results from parameter spaces of ideal oscillators, these Arnold tongue attributes 
are a consequence of the non-ideal character of this oscillator. 

This paper is organized as follows. In Section 2 we present the mathematical descrip¬ 
tion of the non-ideal Duffing oscillator. In Section 3, we show our numerical analysis. 
The last section contains our main conclusions. 


2. Non-ideal Duffing oscillator 


Several mechanical systems can be described by the Duffing equation. Tusset and 
Balthazar ^32] studied ideal and non-ideal Duffing oscillator with chaotic behavior. They 
suppressed the chaotic oscillations through the application of two control signals. In this 
work, we consider a non-ideal system consisting of a mass, spring and pendulum. Figure 
Q] shows a schematic model of the non-ideal oscillator 31], that is composed of a cart 
(mass M), with a pendulum (mass m and length r), connected to a fixed frame by a 
nonlinear spring and a dash-pot. We denote by X the displacement of the cart and by 
t p the angular displacement of the pendulum. 
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Figure 1: Schematic model of the non-ideal oscillator. 


The equations of motion, obtained by using Lagrangian approach, for both 


the cart 


and the pendulum are given by: 


, ^,d 2 X dX , v 

(m + M) —|- ci— 7 -— k\X + JV 2 X 

dt- dt 

od 2 tp dip 

mr -tit- + C 9 ——I- mgr sm p 

dt z dt 


dp‘ 


d 2 p 


mr lift smip-—cosp 


= E — mr 


d 2 X 


dt 2 


cos p, 


( 1 ) 

( 2 ) 


where E is a constant source of energy. According to Equation (]T|). for k\ < 0, the Duffing 
oscillator can be interpreted as a forced oscillator with a spring whose restoring force is 
F = k\X — ^X 3 . Whereas, for k\ > 0, the Duffing oscillator describes the dynamics of 
a point mass in a double well potential, such as a deflection structure building model. 

Considering x = X/r and r = u)\t (wi = m+M )» equations of motion are 
rewritten in the following form: 


x + f3\x — x + "/x 3 = s [ip 2 sin <p — (p cos ip) , 

P + hp + ff 2 sin</? = a — xcosp. 


( 3 ) 

( 4 ) 


for /3i = 




ry = 1=2. '2 

I ~ fcl ' > 


£ = 


m 


m+M : 


h = 


n = St ( w 2 = Vg/r), and 


a ~ Lb (source of energy). 

These equations of motion correspond to a simplified mathematical model for oscil¬ 
lator with a limited power supply. In this case, the source of energy is given by a DC 
motor and the parameter a is associated with its input voltage. 


3. Arnold tongues and Shrimps 

Many systems exhibit behaviour that can be studied by means of two-parameter 


• sv i 

analysis 33,134]. Baesens and collaborators [35| studied two-parameter families of torus 
maps that involve change of mode-locking type. They analysed codimension-1 and -2 
bifurcations for flows on the torus, and found a large variety of bifurcation diagrams, 
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Figure 2: (Color online) (a) Bifurcation diagram showing coexisting attractors for x in terms of a with 
e = 0.09. (b) Magnification of bifurcation diagram for the attractors plotted in black. 


some of them with Arnold tongues. Their numerical studies were focussed on resonance 
regions revealing a rich collection of codimension-one and -two bifurcations. 

We analyze, in this Section, self-similar structures and Arnold tongues in the param¬ 
eter space of the non-ideal Duffing oscillator were performed by using the fourth-order 
Runge-Kutta method with a fixed step. The control parameters were fixed at = 0.05, 
P 2 = 1.5, 7 = 0.1, and D = 1.0. We consider for dynamic investigations the variations of 
parameters e (the ratio of the masses) and a (input voltage of the DC motor). 

First of all, we use a bifurcation diagram, as shown in Figures^ a) and (b) for e = 0.09, 
to verify possible solutions generated by the oscillator. This diagram is constructed 
varying the control parameter a. For each value of the parameter, we plot the local 
maximum values of the dynamical variable x neglecting the transients. As can be seen in 
Figures[D(a) and (b), the bifurcation diagram is composed of periodic windows associated 
with period-adding sequences. See as an example the three main periodic windows in 
(b), with periodic attractor of periods 5, 4, and 3. Then, as a is increased the period 
decreases by 1 . 

In addition, figure [21(a) exhibits the coexistence of two attractors each one plotted 
with a different color (black and blue). In mechanical systems, the coexistence of attrac¬ 
tors is quite common non-linear phenomenon observed. For example, the coexistence of 
a large number of periodic attractors in a mechanical system was observed by Feudel and 
collaborators [36j]. They studied the kicked double rotor system, and verified the possi¬ 
bility of the system to be stabilized by means of a small perturbation. In an experimental 
nonlinear pendulum, it was also observed two coexisting attractors 37j. Multiple attrac- 
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tors may be found in many nonlinear dynamical systems, for instance, driven damped 
pendulum (39} and spring-pendulum system [40]. 



Figure 3: (Color online) Parameter plane diagram for e versus a with different initial conditions. Periodic 
solutions are plotted in blue, green and black scale ranges, quasi-periodic in red, and chaotic in yellow 
and white. 

In order to better characterize the dynamics of the oscillator and to examine the 
structures related to the periodic windows, we construct diagrams of two-dimensional 
parameter space by using the Lyapunov exponents. To evaluate these exponents, we use 
the algorithm proposed by Wolf and coworkers [38^ . One positive Lyapunov exponent 
(LLE) indicates a chaotic attractor, all negative exponents a periodic, and two null 
exponents a quasi-periodic or a bifurcation point. 

In Figures [3])a) and (b), we plot the parameter plane diagrams, for e versus a, using 
a grid of 800x800 cells. Periodic solutions are plotted in blue, green and black scale 
ranges, quasi-periodic in red (bifurcation points in red), and chaotic in yellow and white, 
where the colors, corresponding to range of Lyapunov exponents values, are introduced to 
emphasize the structures details. To make evident the coexistence of multiple attractors, 
these diagrams were constructed considering different initial conditions. From an applied 
perspective, system with coexistence of attractor in noisy environments lead to basin 
hopping [4l|, |42|, the alternate switching among different attractors. 

In Figure^a), we provide a magnification of rectangular area of Figure[3])a) revealing 
many periodic structures (in black) knowing as Arnold tongues f6] that correspond to 
phase locking, i.e., periodic orbits with the same frequency. Surprisingly, the tongues 
origins appear for low value of e for a given a, accumulate in a starting point, namely, 
the tongues distributions appear highly organized with a codimension-2 bifurcation as 
a point of accumulations. In this case, the point of accumulations corresponds to a 
saddle-node Hopf bifurcation [ 35 }. 

Moreover, we can observe in Figures QJb) and (c) small periodic structures (blue, 
green, and black), called shrimps Q, embedded in parameter regions with chaotic regimes 
(yellow and white). The shrimps are composed of the central body bordered by a saddle- 
node and a flip bifurcations. As can be noted in Figure[4])d), the shrimp-shaped windows 
present an interesting feature of self-similar properties. In other words, we can verify in 
successive magnifications (not shown here) the repetition of these types of windows for 
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Figure 4: (Color online) (a) Magnification of the parameter plane diagram shown in Fig. [3ji. (b)-(d) 
Successive magnifications. Periodic solutions are plotted in blue, green and black scale ranges, quasi- 
periodic in red, and chaotic in yellow and white. 


different length scales. 

In the end, we provide in Figure[5ja) a magnification for a rectangular area involving 
the accumulation point of the Arnold tongues. The inferior part of tongues present 
unexpected shape being quite similar to the superior part. In addition, embedded in 
chaotic regions an uncountable number of shrimps as shown in Figure [5](b). 

In addition, the periodic windows in the parameter spaces can be understood as areas 
delimitaded by the codinrension-2 saddle-node Hopf bifurcation curves. On the other 
hand, the Arnold tongues appear whenever the ratio between the driven and natural 
frequencies is a rational number, and they accumulate in a codinrension-2 point. 

4. Conclusions 

We have investigated the parameter space of a non-ideal Duffing oscillator, namely, 
the mass-spring-pendulum system. Initially, we have identified coexisting attractors with 
period-adding cascades. In this case, it is possible to verify periodicities of the windows 
by using the Fibonacci rule [f>]. This rule characterizes an integer sequence, in that the 
sequence is the sum of two previous ones. Moreover, this sequence is related with Golden 
ration. 
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Figure 5: (Color online) (a) Magnification of parameter plane diagram shown in Fig. |3ji. (b) Magnifica¬ 
tion of (a). 


When the source of energy was included in the oscillator we were able to observe 
parameter regions identified as Arnold tongues corresponding to mode locked and pe¬ 
riodic motion with a common frequency. The mode locked occurs when the combined 
motion presented in the mass-spring-pendulum driven by a DC motor becomes periodic. 
We have verified that locked and unlocked regions were interwoven in parameter space. 
Our results showed that the tongues origins accumulate in a low value of e for a given 
a. Furthermore, the tongues are organized from a codimension-2 bifurcation as a point 
of accumulations. We also observed shrimp-shaped structures immersed in parameter 
regions with chaotic regimes, and with highly organized distributions. We have found 
shrimps in small range of the parameter space. 

Moreover, it is important to emphasize that characterization of the attractors in 
parameter space of applied systems is useful to choose robust periodic orbits and also to 
evaluate the attractor changes in case of controlling chaotic oscillations. 
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